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Abstract.  We  measure  high- resolution  Feshbach  resonance  spectra  for  ultracold  ce¬ 
sium  atoms  colliding  in  different  hyperfine  and  magnetic  sublevels.  More  than  25 
Feshbach  resonances  are  observed  for  magnetic  fields  below  230  G  in  the  elastic  and 
inelastic  ground  state  collision  cross  sections,  as  well  as  in  the  cross-section  for  light- 
assisted  collisions.  From  these  spectra  a  consistent  set  of  ground  state  molecular  in¬ 
teraction  parameters  for  cesium  is  extracted,  including  singlet  and  triplet  scattering 
lengths  of  As  =  (280±10)ao  and  At  =  (2400±100)ao,  a  van  der  Waals  coefficient  Ce 
=  (6890±35)a.u.,  as  well  as  the  strength  of  the  indirect  spin-spin  coupling.  This  set  of 
parameters  allows  for  the  first  time  a  complete  characterization  of  cesium’s  ultracold 
collision  properties. 


Ultracold  cesium  collisions  have  been  the  subject  of  a  large  number  of  experi¬ 
mental  [1-8]  and  theoretical  [9-12]  studies  in  recent  years.  Aside  from  its  obvious 
importance  as  the  source  of  our  primary  time  and  frequency  standard,  recent  ex¬ 
periments  have  shown  that  ground  state  cesium  exhibits  anomalous  collision  prop¬ 
erties,  such  as  a  zero-energy  resonance  [2]  and  extremely  large  dipolar  relaxation 
rates  [3,4] .  These  effects  have  important  implications  for  the  stability  of  laser-cooled 
cesium  atomic  clocks  [5],  precision  measurements  such  as  searches  for  T-violating 
effects  beyond  the  standard  model  in  atoms  [13],  as  well  as  attempts  to  produce  a 
cesium  Bose- Einstein  condensate  [4]. 

In  this  paper  we  present  high-resolution  studies  of  these  cold  collision  interactions 
using  a  technique  we  label  Feshbach  spectroscopy.  This  method  exploits  the  fact 
that  the  threshold  collision  properties  of  a  binary  system  are  almost  completely 
determined  by  the  presence  of  molecular  bound  or  virtual  states  near  the  threshold 
for  dissociation.  We  probe  these  very  weakly  bound  levels  by  tuning  them  into  de- 
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generacy  with  the  scattering  continuum  of  a  lower  lying  internal  state  of  the  atom 
pair  using  the  relative  Zeeman  shift  induced  by  an  external  magnetic  held.  When  a 
bound  level  is  resonant  with  colliding  atom  pairs,  the  elastic  and  inelastic  collision 
cross-sections  exhibit  a  resonance  behavior  hrst  predicted  by  H.  Feshbach  in  the 
context  of  nuclear  scattering  [14] ,  and  more  recently  by  Verhaar  and  co-workers  for 
cold  atomic  collisions  [10] .  These  so-called  Feshbach  resonances  are  ubiquitous,  and 
have  been  either  predicted  [10,11,15]  or  observed  [6-8,16,17]  to  occur  in  virtually 
all  of  the  alkali  atomic  species.  Since  Feshbach  resonance  spectroscopy  involves 
only  atoms  in  stable  ground  states,  it  has  inherently  high  resolution,  allowing  us 
to  precisely  chart  the  positions  of  weakly  bound  molecular  states,  and  from  these 
measurements  extract  for  the  hrst  time  a  self-consistent  set  of  ground  state  molecu¬ 
lar  interaction  parameters  for  cesium  [8,12].  Possible  implications  for  future  atomic 
clocks  and  cesium  Bose-Einstein  condensation  are  discussed. 

The  apparatus  used  for  these  measurements  has  been  described  previously  [6-8] . 
Samples  of  up  to  2  x  10®  cesium  atoms  are  loaded  into  a  far-detuned  dipole  trap 
from  a  vapor-cell  MOT  using  3D  degenerate  Raman  sideband  cooling  in  a  3D  opti¬ 
cal  lattice  [18]  which  is  superimposed  on  the  overlapping  MOT  and  dipole  trapping 
volume.  Temperatures  in  the  range  of  a  few  pK  and  densities  above  10^®cm^®  are 
produced  in  the  dipole  trap  which  is  formed  by  a  20W  single  frequency  Nd:YAG 
laser  retrorehected  upon  itself  to  form  a  ID  optical  lattice,  which  is  sufficiently  far 
detuned  from  resonance  to  make  the  spontaneous  scattering  rate  of  trap  photons 
negligible.  The  lifetime  of  the  atoms  in  this  trap  is  background-gas  pressure  limited 
to  about  2  seconds,  providing  ample  time  for  sensitive  collisional  measurements. 
The  number  of  atoms  in  the  trap  is  measured  by  collecting  resonance  huorescence 
onto  a  calibrated  photodiode,  the  temperature  of  the  sample  is  measured  using  a 
time-of-hight  technique,  and  the  trap  vibration  frequencies  using  parametric  exci¬ 
tation  [6-8]. 

We  observe  collisional  processes  in  atomic  samples  polarized  into  several  diherent 
magnetic  sublevels  of  the  cesium  ground  state.  In  all  cases  the  atomic  polarization 
is  optimized  with  care,  and  polarization  purities  as  high  as  97%  are  achieved.  First, 
we  consider  |F=3,m=3),  into  which  the  atoms  are  already  polarized  after  Raman- 
sideband  cooling.  This  is  the  lowest  energy  state  in  a  magnetic  held,  and  therefore 
it  is  stable  against  binary  inelastic  processes.  We  observe  Feshbach  resonances 
in  the  elastic  collisions  between  these  atoms  by  loading  them  into  a  shallow  trap 
whose  depth  is  comparable  to  their  temperature,  and  monitoring  the  subsequent 
evaporative  loss  from  the  trap  as  a  function  of  magnetic  held.  A  large  elastic 
cross-section  results  in  a  fast  loss  of  atoms  from  the  trap,  while  a  small  elastic 
cross-section  produces  a  slow  loss.  Fig.  1(a)  shows  this  evaporative  loss  over  a 
hxed  time  interval  as  a  function  of  magnetic  held.  Although  this  quantity  is  not 
a  linear  function  of  the  elastic  cross-section,  it  is  still  a  monotonic  function  and 
provides  higher  signal-to-noise  for  the  location  of  resonance  features  in  the  cross- 
section  than  previous  methods  using  temperature  measurements  [6]. 

To  prepare  atoms  in  the  |3,-3)  state,  which  is  of  most  interest  for  magnetic 
trapping  and  evaporative  cooling  [3,4],  we  quickly  reverse  the  magnetic  bias  held 


magnetic  field  [G] 


FIGURE  1.  T(a)  Elastic  collision-induced  loss  rate  from  a  shallow  trap  vs.  magnetic  field  for 
the  |3,3)  state,  which  is  stable  against  inelastic  binary  collisions,  (b)  Theoretical  curve  obtained 
by  fitting  to  the  resonance  positions  shown  in  Figs.  1(a),  2,  and  3.  The  positions  of  the  minima 
exactly  match  those  in  Fig.  1(a). 


along  which  the  atoms  are  polarized  during  sideband  cooling,  projecting  a  large 
fraction  of  the  atoms  immediately  into  the  desired  state.  This  is  followed  by  a  short 
optical  pumping  pulse.  As  was  already  observed  in  Ref.  [3],  this  state  decays  quickly 
via  dipolar  relaxation  (arising  from  the  magnetic  dipole-dipole  interaction  between 
the  spins  of  the  valence  electrons)  to  other  sublevels  within  the  F=3  manifold.  Note 
that  this  state  is  stable  against  inelastic  spin-exchange  collisions,  since  in  order  to 
decay  it  must  change  its  total  angular  momentum  along  the  magnetic  held,  which 
cannot  happen  via  the  spherically  symmetric  ground-state  exchange  interaction. 
The  magnetic  potential  energy  released  in  the  decay  is  typically  much  larger  than 
the  depth  of  our  dipole  trap,  so  that  the  inelastic  collision  rate  coefficient  can  be 
easily  extracted  from  the  observed  time-dependent  trap  loss  rate.  This  quantity  is 
shown  in  Fig.  2  vs.  magnetic  held. 

Note  the  large  number  of  Feshbach  resonances,  in  comparison  to  the  other  alkalis 
where  at  most  one  or  two  have  been  predicted  or  observed  in  a  magnetic  held  range 
of  this  size  [15-17].  This  is  due  to  the  much  stronger  inhuence  of  the  so-called 
indirect  spin-spin  coupling,  where  the  ground  state  valence  electron  spins  are  mag¬ 
netically  coupled  via  the  spin-orbit  interaction  to  distant  excited  electronic  states 
of  the  dimer  having  orbital  angular  momentum  [19].  The  result  is  a  coupling  of 
the  total  electron  spin  to  the  molecular  axis  which  breaks  the  spherical  symme- 


FIGURE  2.  Dipolar  loss  rate  coefficient  vs.  magnetic  field  for  atoms  in  the  |3,-3)  state.  The 
solid  line  is  the  theoretical  curve  obtained  by  fitting  to  the  resonance  positions  shown  in  Figs. 
1(a),  2,  and  3. 


try  of  the  ground  state  molecular  interaction  in  the  same  way  as  the  well-known 
(’’direct”)  magnetic  dipole-dipole  interaction  [19],  and  therefore  contributes  to  the 
cross-section  for  dipolar  relaxation.  For  the  lighter  alkalis  this  is  an  extremely  small 
effect,  and  the  already  very  weak  direct  contribution  dominates  these  processes, 
however  for  Cs2  the  indirect  contribution  is  anomalously  large  and  overwhelms  the 
direct  coupling  [19,11,12].  This  has  two  important  consequences;  hrst,  the  dipolar 
relaxation  cross-sections  are  orders  of  magnitude  larger  than  for  the  lighter  alkalis 
[3,4].  Second,  free  Cs  atom  pairs  in  S-wave  scattering  states  can  couple  much  more 
strongly  to  rotating  molecular  states  of  the  dimer,  and  this  produces  the  many  nar¬ 
row  Feshbach  resonances  observed  in  these  relaxation  rates.  The  positions  of  these 
resonances  provide  a  direct  measurement  of  the  strength  of  the  indirect  spin-spin 
coupling,  which  has  been  an  important  source  of  uncertainty  in  the  Cs-Cs  interac¬ 
tion,  and  has  until  now  helped  to  prevent  a  consistent  picture  of  these  interactions 
from  emerging  from  the  available  experimental  and  theoretical  information. 

Finally,  we  observe  collisions  between  atoms  polarized  in  the  |4,-4)  state.  The 
observed  rate  coefficient  vs.  magnetic  held  is  shown  in  Fig.  3.  These  samples  are 
prepared  using  a  fast  held  reversal  as  described  above,  again  followed  by  optical 
pumping.  Loss  from  this  state  occurs  mainly  via  hyperhne-changing  collisions, 
where  one  or  both  of  the  atoms  decays  through  dipolar  relaxation  to  the  lower 
F=3  manifold.  Resonances  in  this  state  provide  an  important  constraint  on  the 
interaction  parameters  which  is  complementary  to  those  already  described,  since 
collisions  in  this  state  involve  only  the  triplet  potential,  whereas  atoms  in  the 


FIGURE  3.  Inelastic  loss  rate  coefficient  vs.  magnetic  field  for  atoms  in  the  |4,-4)  state.  The 
solid  line  is  the  theoretical  curve  obtained  by  fitting  to  the  resonance  positions  shown  in  Figs. 
1(a),  2,  and  3. 


lower  F=3  manifold  experience  a  mixture  of  the  singlet  and  triplet  potentials.  In 
addition,  the  bound  states  responsible  for  Feshbach  resonances  in  collisions  of  |4,-4) 
atoms  must  have  binding  energies  less  than  the  Zeeman  energy  of  the  atom  pair, 
rather  than  the  hyperhne  energy  as  is  the  case  for  atoms  colliding  in  the  lower  F=3 
manifold.  These  extremely  weakly  bound  states  are  therefore  more  sensitive  to  the 
long-range  part  of  the  potential  [19]. 

We  use  a  coupled-channel  collision  model  similar  to  that  described  in  Refs  [11,19] 
to  ht  the  observed  resonance  positions  shown  in  Figs.  1(a),  2,  3  [12].  The  input  to 
this  model  are  the  following  four  htting  parameters:  the  singlet  and  triplet  scat¬ 
tering  lengths  Ag  and  At,  the  van  der  Waals  coefficient  Cq,  and  a  dimensionless 
scaling  factor  Sc  which  controls  the  strength  of  the  indirect  spin-spin  interaction 
[11,12].  The  hts  are  shown  as  a  solid  line  in  Figure  1(b)  and  in  Figs.  2-3  superim¬ 
posed  on  the  data,  and  they  reproduce  the  positions  of  all  the  observed  resonances. 
The  resulting  parameter  values  are  Ag  =  (280  -t— 10)ao,  At  =  (2400+-100)ao,  Cq 
=  (6890-1— 35)a.u.,  and  Sc  =  3.2-(— 0.5.  Notice  that  both  the  singlet  and  triplet 
scattering  lengths  are  positive  and  very  large;  these  arise  from  extremely  weakly 
bound  states  in  the  singlet  and  triplet  potentials,  the  latter  of  which  is  the  source 
of  the  zero  energy  resonance  observed  in  the  pure  triplet  state  [2] . 

In  addition  to  the  three  spectra  shown  in  Figs.  1(a),  2,  and  3,  which  were  taken 
with  spin-polarized  samples  and  were  used  to  extract  values  for  the  interaction  pa¬ 
rameters,  we  have  also  observed  many  additional  Feshbach  resonances  in  collisions 
between  atoms  in  different  magnetic  sublevels,  where  P-wave  collisions  are  possi- 
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FIGURE  4.  Inelastic  loss  rate  coefficient  vs.  magnetic  field  for  a  mixed  sample  with  70%  |3,-3), 
25%  |3,-2),  and  5%  |3,-1).  Note  that  many  of  the  additional  features  are  due  to  P-wave  collisions 
between  atoms  in  different  sublevels.  The  solid  line  is  the  theoretical  curve  obtained  by  fitting  to 
the  resonance  positions  shown  in  Figs.  1(a),  2,  and  3. 

ble.  A  detailed  presentation  of  these  spectra  and  resonance  positions  is  beyond  the 
scope  of  this  paper,  and  will  be  presented  elsewhere.  However,  we  give  an  exam¬ 
ple  in  Fig  4.  This  hgure  shows  the  observed  inelastic  loss  for  a  sample  of  atoms 
in  a  mixture  of  the  |3,-3),  |3,-2),  and  |3,-1)  states,  which  was  obtained  simply  by 
imperfectly  polarizing  the  gas.  All  additional  resonances  are  correctly  reproduced 
by  the  calculations,  whose  only  input  is  the  four  parameters  extracted  from  Figs. 
1(a),  2,  and  3. 

We  now  call  the  reader’s  attention  to  the  small  resonances  predicted  in  the  elastic 
cross  section  for  |3,3)  atom  collisions,  which  are  magnihed  in  the  inset  of  Fig.  1, 
and  are  clearly  too  weak  to  be  seen  in  the  evaporative  loss.  In  order  to  observe 
these  features,  we  use  optical  excitation  to  probe  the  short-range  molecular  forma¬ 
tion  characteristic  of  Feshbach  resonance  [17].  A  free  atom  pair  within  an  optical 
wavelength  of  each  other  experiences  a  very  strong,  spatially  dependent  shift  of 
the  electronic  transition  frequency  due  to  the  resonant-dipole  interaction  in  the 
excited  state.  A  light  held  at  a  certain  detuning  will  thus  be  resonant  only  at  a 
particular  interatomic  distance,  called  the  Condon  point  Rc,  and  the  excitation 
probability  is  proportional  to  the  probability  density  with  which  atoms  are  found 
at  that  particular  separation  [20].  Excited  pairs  can  undergo  a  so-called  radiative 
collision,  where  the  tremendously  strong  dipole  forces  in  the  excited  state  acceler¬ 
ate  the  atoms  to  a  sufficiently  high  kinetic  energy  before  they  decay  that  they  are 
expelled  from  the  trap  [20,7].  By  monitoring  the  light-induced  loss  from  our  trap 
at  a  suitable  detuning  (in  this  case  6  nm  above  resonance),  we  can  directly  probe 
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FIGURE  5.  (a)  Radiative  collision  loss  rate  vs.  magnetic  field  for  atoms  in  the  |3,3)  state,  (b) 
Calculated  scattering  length  showing  the  predicted  Feshbach  resonances.  Note  that  even 
the  very  weak  features  are  clearly  seen  in  the  loss  rate.  Several  of  the  observed  resonances  are 
not  predicted,  and  may  be  due  to  higher-oder  scattering  processes. 


atomic  pair  correlation  function  at  a  small  interatomic  distance  characteristic  of 
molecular  formation,  and  thus  sensitively  detect  the  resonant  enhancement  of  this 
quantity  which  occurs  at  a  Feshbach  resonance.  Figure  5(a)  shows  the  observed 
radiative  collision  loss  for  atoms  in  the  |3,3)  state  as  a  function  of  magnetic  held, 
and  Figure  5(b)  shows  the  predicted  scattering  length  of  this  state  for  comparison 
of  the  resonance  positions.  Both  of  the  weak  features  predicted  by  the  calculations 
appear  clearly  in  the  radiative  loss  rate,  as  well  as  the  larger  resonance  at  48  G 
already  observed  in  the  evaporative  loss.  Note  also  that  several  additional  features 
are  also  present  which  were  not  predicted  by  the  theory.  We  are  currently  investi¬ 
gating  the  possibility  that  these  resonances  result  from  higher  order  processes  where 
the  incoming  S-wave  couples  to  rotating  molecular  states  of  even  higher  angular 
momentum  (e.g.  1=4,6  etc...),  and  the  agreement  appears  to  be  promising.  This 
will  be  the  subject  of  a  future  publication. 

With  a  consistent  set  of  ultracold  molecular  interaction  parameters,  we  can  now 
calculate  the  rates  for  arbitrary  cold  collision  processes.  Of  particular  interest  is 
the  cold-collision  induced  frequency  shift  of  the  so-called  ’’clock  transition”  between 
|3,0)  and  |4,0),  which  is  likely  to  be  the  most  important  limitation  on  the  short-term 
stability  of  the  next  generation  of  laser-cooled  cesium  atomic  fountain  clocks  [5,21]. 


FIGURE  6.  (a)  Cold-collision  induced  frequency  shift  of  the  |3,0)  to  |4,0)  transition  calculated 
for  an  equal  mixture  of  these  two  states,  as  a  function  of  collision  energy.  The  existing  experimental 
points  are  shown,  and  agree  with  the  calculation.  Note  that  the  shift  goes  through  zero  at  about 
90  nK. 


Shown  in  Fig.  6  is  the  calcnlated  freqnency  shift  per  nnit  density,  as  a  fnnction  of  the 
atomic  kinetic  energy  in  temperatnre  nnits.  Also  shown  are  several  experimental 
points  from  Ref.  [5],  which  are  in  good  agreement  with  the  calcnlation.  Notice 
that  the  cnrve  passes  throngh  zero  and  changes  sign  at  a  temperatnre  of  abont 
90  nK.  If  a  large  enongh  sample  of  atoms  conld  be  prodnced  in  a  narrow  energy 
range  aronnd  this  point,  it  might  be  possible  to  dramatically  rednce  the  effect  of 
this  density-indnced  shift  [21].  3D  Raman-sideband  cooling  [18]  has  already  been 
nsed  to  prodnce  an  atomic  fonntain  with  10®  atoms  at  temperatnres  as  low  as  150 
nK  [22],  and  wonld  therefore  be  a  very  promising  techniqne  for  this  pnrpose. 

As  a  second  and  hnal  example  of  what  can  be  calcnlated  nsing  the  interaction 
parameters,  we  now  consider  the  impact  this  new  information  has  on  attempts  to 
prodnce  a  cesinm  Bose-Einstein  condensate.  Shown  in  Fig.  7  are  the  calcnlated 
elastic  and  inelastic  collision  rate  coefficients  for  atoms  in  the  |3,-3)  state,  where 
experimental  efforts  at  BEC  in  a  magnetic  trap  have  been  focnsed  [3,4].  The  ratio 
of  these  two  qnantities  R  is  a  hgnre  of  merit  for  forced  evaporative  cooling  [23]. 
Althongh  R  takes  on  valnes  of  100  or  larger  away  from  the  resonances,  there  are 
additional  considerations  that  are  important.  First,  the  sign  of  the  scattering  length 
mnst  be  positive,  at  least  dnring  the  hnal  evaporation  stages  where  T<Tc,  to  have 


FIGURE  7.  For  atoms  in  the  |3,-3)  state,  the  calculated  elastic  collision  rate  coefficient  is  shown 
in  (a),  the  inelastic  collision  rate  coefficient  in  (b),  and  the  scattering  length  in  (c).  The  collision 
rate  coefficients  are  orders  of  magnitude  larger  than  typically  seen  in  the  lighter  alkalis;  also,  the 
scattering  length  for  |3,-3)  is  large  and  negative  in  the  low- field  region  where  evaporation  was 
performed  in  Ref  [4], 


a  large,  stable  condensate.  This  eliminates  the  region  below  about  lOG,  where  the 
previous  experiments  were  performed  [3,4],  Second,  the  absolute  magnitudes  of 
these  rate  coefficients  are  enormous.  This  is  important  because  the  trap  vibration 
period  and  trap  lifetime  are  hxed  timescales  against  which  elastic  and  inelastic 
rates  must  be  compared  absolutely,  rather  than  just  to  each  other.  Inelastic  loss 
rates  must  be  kept  small  enough  so  that  one  has  sufficient  time  to  observe  the 
condensate,  and  elastic  collision  rates  larger  than  the  trap  vibration  frequency 
cause  a  crossover  into  hydrodynamic  behavior.  The  latter  may  be  a  serious  problem 
for  evaporation,  either  by  limiting  the  thermalization  rate  [6]  or  by  reducing  the 
evaporation  efficiency  or  dimensionality.  These  effects  must  be  characterized  before 
an  intelligent  choice  of  parameters  can  be  made  based  upon  our  measurements. 

In  conclusion,  we  have  presented  high-resolution  spectra  of  elastic  and  inelastic 


collision  rates  in  an  nltracold,  trapped  cesinm  gas.  From  these  spectra  we  have 
for  the  hrst  time  extracted  a  self-consistent  set  of  molecular  interaction  parameters 
which  reproduce  the  positions  of  over  30  observed  Feshbach  resonances.  These  pa¬ 
rameters  can  be  used  to  accurately  calculate  important  collisional  effects  such  as 
the  density-dependent  frequency  shift  in  cesium  atomic  clocks,  as  well  as  collision 
rates  and  scattering  lengths  relevant  for  attempts  at  Bose-Einstein  condensation. 
In  addition,  Feshbach  spectroscopy  may  be  more  widely  useful  as  a  general  spec¬ 
troscopic  tool  for  studying  cold  collisions,  especially  when  optical  excitation  is  used 
as  a  sensitive  probe  of  the  molecular  formation. 
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